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INFLUENCE  OF  RING  STIFFENERS  ON  INSTABILITY  OF 
ORTHOTROPIC  CYLINDERS  IN  AXIAL  COMPRESSION* 

By  David  L.  Block 
Langley  Research  Center 

SUMMARY 

f .  >  __ ) 

1  Calculations  are  presented /from  an  analytical  investigation  on  the  influence 
of  ring  stiffeners  on  the  insta^lity  modes  of  orthotropic  cylinders  subject  to 
compressive  or  bending  loads.  pThe  analysis  is  performed  by  employing  small- 
deflection  theory  and  by  modifying  the  equilibrium  equati^^^  include  the 
effects  of  discretq>.ri|^  stiffeners  characterized  by  a  bending  stiffness  that 
restrains  radial  dejbfeation  of  the  shell  calculations  indicate  that  the 

ring  bending  stiffness  necessary  to  cau'^  panel  instability  can  be  adequately 
determined  by  use  of  an  analysis  which  does  not  include  the  discreteness  of  the 
rings.  Comparison  of  the  results  of  the  calculations  with  an  empirical  ring- 
design  criterion  in  common  use  indicates  that  the  eiiq)irical  formula  can  be  either 
very  conservative  or  very  nonconservative  depending  on  the  cylinder  geometr^^  |  ^ 


INTRODUCTION 

encountered  in  the  design  of  axially  stiffened,  axially  com¬ 
pressed  cylinders  is  the  determination  of  the  size  of  the  circximferential  stiff¬ 
ening  elements  or  rings  required  to  prevent  general  instability  failure  of  "the 
cylinders.  Common  practice  is  to  determine  the  size  of  rings  by  the  empirica 
formula  lof  reference  1, [which  gives  the  required  ring  bending  stiffness  to  force 
cylinder  failure  to  occur  between  rings_^  The  formula  of  reference  1  is  base 
on  tests  of  small  cylinders  with  relatively  few  stiffening  elements.  Such 
cylinders  are  not  very  representative  of  the  cylinders  used  in  contemporary  air¬ 
craft  and  launch  vehicles.  The  use  of  this  formula  in  contemporary  design 
therefore  causes  some  concern. 

fAn  analytical  study  of  the  ring  stiffness  required  to  force  cylinder  fail¬ 
ure  to  occur  between  rings  is  made  in  the  present  study^^  For  this  purpose  an 
analysis  incorporating  a  Donnell  type  theory  similar  to  the  one  employed  in 
references  2  and  3  was  used,  p^e  mathematical  treatment  of  the  rings  is  approx¬ 
imate  and  the  ring  stiffnesses  predicted  by  the  calculations  madejh-e^ln  cannot 
be  considered  conclusive.  However,  the  calculations  do  Indicate 
acteristics  which  aid  in  the  design  of  torsionally  weak  rings  of  axiallyj 

- *The  information  presented  herein  is  based  in  part  upon  a  thesis  offered  in 

partial  fulfillment  of  the  requirements  for  the  degree  of  Master  of  Science  in 
Engineering  Mechanics,  Virginia  Polytechnic  Institute,  Blacksburg,  Virginia, 

June  1964- 


rstiffened  cylinders  and  which  give  a  ring  design  criterion  that  has  a  better 
theoretical  basis  than  the  one  i  (ref.  l)yiibw  in  general  use^ 

^ —  i<i  p- 


SYMBOLS 
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a 


Dx,Dy 


Dxy 

d 

El* 


Gxy 


elements  of  a  determinant 

coefficient  of  deflection  function 

transverse  shear  stiffnesses  of  cylinder  wall  in  longitudinal  and 
circumferential  directions,  respectively 

bending  stiffness  of  cylinder  wall  in  longitudinal  and  circumferential 
directions,  respectively 

twisting  stiffness  of  cylinder  wall 

stiffener  spacing,  ^  ^ 

Young's  modxilus  of  elasticity  for  ring  stiffeners 

extenslonal  stiffness  of  cylinder  wall  in  longitudinal  and  circumfer¬ 
ential  directions,  respectively 

In-plane  shear  stiffness  of  cylinder  wall 


H  = 


Dx 


Ip  moment  of  inertia  of  ring  stiffeners 

i,  integers 


kx  critical-axial- compressive- stress  coefficient^ 

L  length  of  cylinder 


- ^  ^  +  2Dxy  +  — + 


'E 


%  Sx^ 


*^xy\  ^  ^xy 


Dy 

1  -  ^^x^V 


inverse  of  Lj; 


defined  by 


2 


Mp  = 


1  -  l^xMy 


+  25^32  + 


dTTi 


bc{-*-  “  ^^x^‘y 


p^z2 


^  fe  -  ^  I 


critical  resiiltant  normal  force  in  axial  direction 

n  number  of  rings 

R  radius  of  cylinder  to  midplane 

s  minor-determinant  nvimber  whose  value  is  one  through  number  of  rings, 

1,  2,  3^  •  •  •  n 

w  displacement  in  radial  direction  of  middle  surface  of  cylinder 

X  longitudinal  coordinate  of  cylinder 

y  circumferential  coordinate  of  cylinder 


cylinder  curvature  parameter. 


R  yDx 


ratio  of  cylinder  length  to  circumferential  buckle  length,  L/A 

ratio  of  bending  stiffness  of  ring  to  bending  stiffness  of  cylinder 
wall  in  circumferential  direction,  EpI^^yd 

ratio  of  bending  stiffness  of  ring  to  bending  stiffness  of  cylinder 
wadi  in  circumferential  direction  at  which  buckling  in  the  panel 
and  general  instability  modes  coincide  ® 

/oo 

f(x)  8(x  -  id)dx  =  f(id), 

.00 

where  5(x  -  Id)  =  0  when  x  ^  id 
Kronecker  delta,  =  0  when  i  /  j  =  1  when  i  =  j 

circumferential  buckle  length 

Poisson’s  ratios  associated  with  bending  of  cylinder  wall 


^^x.>  My 
^lx^My’ 


Poisson’s  ratios  associated  with  extension  of  cylinder  wall 
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RESULTS  MD  DISCUSSION 


The  calculations  and  results  presented  were  made  with  the  use  of  the  sta¬ 
bility  criterion  presented  in  appendix  A.  The  criterion  applies  to  orthotropic 
cylinders  loaded  in  compression  and  is  based  on  conventional  small-deflection 
buckling  theory  (ref.  4);  the  rings  are  treated  as  discrete  members  located  at 
the  neutral  axis  of  the  cylinder;  they  are  characterized  by  a  single  stiffness, 
a  bending  stiffness  that  restrains  radial  deformations  of  the  shell  (refs.  2 
and  5)*  Calculations  vere  made  for  cylinders  covering  a  vide  range  of  ortho¬ 
tropic  stiffnesses  in  order  to  study  the  influence  of  the  discreteness  of  rings 
on  the  instability  modes  of  stiffened  cylinders.  Although  the  calculations  were 
made  only  for  cylinders  loaded  in  compression,  they  vill  also  apply  with  reason¬ 
able  accuracy  to  cylinders  subject  to  bending  loads. 

All  results  of  the  calculations  made  exhibited  similar  characteristics  and, 
therefore,  only  a  typical  calculation  is  presented  in  tables  I  and  II  where 
theoretical  buckling  coefficients  for  a  selected  stiffened  cylinder  are  given 
in  terms  of  the  natural  parameters  of  the  problem.  The  data  of  tables  I  and  II 
are  also  given  in  figures  1  and  2  in  the  form  of  plots  of  the  buckling  coeffi¬ 
cient  against  the  ring  stiffness  parameter.  Separate  curves  are  given  for  dif¬ 
ferent  values  of  the  curvature  parameter  and  for  different  modes  of  buckling. 

The  solid  curves  of  figure  1  denote  buckling  of  the  cylinder  in  the  general 
instability  mode  and  the  short  dashed  lines  denote  buckling  in  the  panel  insta¬ 
bility  mode.  It  vill  be  noted  that  the  curves  denoting  general  Instability  form 
an  envelope  curve  which  is  independent  of  the  discreteness  of  the  rings.  This 
curve  can  be  calculated  by  conventional  orthotropic  analysis  (refs.  5  sind  6), 
which  distributes  the  bending  stiffness  of  the  rings  over  the  panel  length 
between  rings  and  treats  it  as  an  additional  property  of  the  cylinder  wall. 

Note  also  that  the  effect  of  ring  discreteness  does  not  become  important  until 
the  buckling  coefficient  for  panel  instability  has  been  exceeded,  except  for  the 
case  of  a  cylinder  with  a  single  ring  (n  =  l)  and  even  for  this  case  the  effect 
is  negligible.  Hence,  cylinders  can  be  adequately  analyzed  by  theories  which 
do  not  account  for  the  discreteness  of  rings.  Another  interesting  feature  of 
figure  1  is  that  the  general  instability  curve  for  a  given  cylinder  levels  off 
at  a  value  of  the  buckling  coefficient  which  is  less  than  the  buckling  coeffi¬ 
cient  for  panel  instability  of  the  cylinder  with  one  additional  ring.  The 
buckling  mode  entails  deflection  of  the  discrete  rings  at  low- values  of  the  ring 
stiffness  ratio  but  as  the  ring  stiffness  ratio  increases  the  deflection  of  the 
rings  approaches  zero  and  buckling  entails  only  deflection  of  the  cylinder  wall 
between  rings. 

Additional  curves  similar  to  those  of  figure  1  are  given  in  figure  2, 
except  that  the  individual  curves  denoting  buckling  by  general  instability  are 
not  included  when  they  do  not  represent  the  governing  mode  of  buckling.  These 
curves  represent  the  governing  mode  of  buckling  only  for  the  cylinder  with  a 
single  central  ring.  (See  curves  for  n  =  1.)  The  general  features  of  the 
curves  of  figure  2  are  the  same  as  those  of  figure  1  of  reference  2  which  gives 
the  buckling  coefficient  of  ring- stiffened  isotropic  cylinders  in  torsion.  In 
both  investigations  the  buckling  coefficient  increases  with  increasing  ring 


k 


s'tiffii6ss  unliil  the  panel  instability  mode  is  reached^  a  id  the  transition  to 
the  panel  instability  mode  occurs  by  a  change  in  mode,  so  that  little  or  no 
Interaction  between  the  two  modes  of  buckling  is  involved.  This  last  feature 
is  important  because  it  Indicates  that  cylinders  can  be  adequately  designed  with 
an  analysis  such  as  the  one  of  reference  6,  which  does  not  include  the  discrete¬ 
ness  of  the  rings.  Both  instability  modes  of  concern  here  are  adequately  dis¬ 
cussed  in  reference  6,  where  for  the  panel  instability  calculation  the  cylinder 
is  considered  to  have  a  length  equal  to  the  ring  spacing  and  for  the  general 
instability  calculation  the  properties  of  the  ring  are  added  to  those  of  the 
wall  in  determining  stiffness  constants  for  the  stability  equation. 

Other  calculations  similar  to  those  presented  in  figures  1  and  2  and  cov¬ 
ering  a  wide  range  of  cylinder  wall  stiffnesses  were  made  in  this  investigation. 
The  resixlts  of  the  calculations  were  similar  to  those  presented  for  a  particular 
set  of  wall  geometries  and  substantiate  the  conclusions  drawn  therefrom. 

A  coD^iarison  of  the  results  of  the  present  theory  with  Shanley's  criterion 
for  ring  stiffness  (ref.  l)  is  given  in  figure  5  for  the  parti c\aar  set  of  wall 
stiffnesses  used  in  the  calculations  for  figures  1  and  2.  Figure  3  was  con- 

2 

structed  with  the  parameter  ( ^  ■■■  held  constantj  this  procedure  results 

in  a  constant  value  for  the  reinforcement  ratio  when  computed  by  the  pres¬ 

ent  analysis.  Figure  3  indicates  that  Shanley's  criterion  may  be  either  very 
conservative  or  very  nonconservative  compared  with  the  present  analysis  depending 
upon  the  proportions  of  the  cylinder.  Additional  comparisons  similar  to  the  one 
shown  in  figure  3  have  been  made  for  cylinders  with  other  wall  stiffnesses. 
Although  the  relative  position  of  the  two  curves  changes,  the  trend  indicated 
in  figure  3  was  found  to  be  typical.  Shanley's  criterion  was  found  to  be  more 
nonconservative  at  the  larger  values  of  d/R  which  correspond  to  small  values 

of  ^  Such  proportions  result  in  a  small  number  of  circumferential  waves 

when  the  cylinder  buckles.  Hence  it  is  concluded  that  Shanley's  criterion  is 
nonconservative  for  cylinders  which  buckle  into  modes  with  a  small  number  of 
circumferential  buckles. 


CONCnXDING  REMARKS 


Resiats  of  theoretical  calculations  investigating  the  influence  of  discrete 
ring  stiffness  on  the  instability  of  orthotropic  cylinders  in  coii:5)resslon  or 
bending  have  been  presented  and  discussed.  A  ring- stiffened  cylinder  can  be 
strengthened  by  increasing  the  ring  stiffness  up  to  a  limit  at  which  the  cyl¬ 
inder  falls  between  the  rings.  The  limit  or  transition  from  general  to  panel 
instability  is  essentially  independent  of  the  discreteness  effect  of  the  rings 
and  allows  the  cylinder  to  be  adequately  analyzed  and  designed  by  conventional 
orthotropic  analysis  which  distributes  the  properties  of  the  rings.  A 


comparison  of  the  analysis  vith  an  empirical  ring  criterion  of  Shanley’s  shows 
that  Shanley’s  criterion  becomes  nonconservative  for  cylinders  which  buckle 
into  modes  with  a  small  number  of  circumferential  buckles* 


Langley  Research  Center, 

National  Aeronautics  and  Space  Administration, 

Langley  Station,  Hampton,  Va.,  August  IJ,  196^. 
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APPENDIX  A 


THEORETICAL  SOLUTION 


The  equation  of  equilibrium  governing  the  buckling  of  a  cu^ed  orthotropic 
plate  with  infinite  transverse  shear  stiffnesses  =  Dq^  =  coj  and  subject 

to  in-plane  axial  stresses  only  is  (ref.  4) 


LdW 


'xy 


Le 


-1  ^ 


—  6(x  -  id)  =  0 


(Al) 


The  fourth  term  of  equation  (Al)  is  added  to  represent  the  radial  restoring 
force  due  to  the  bending  stiffness  of  equally  spaced  rings.  (See  refs.  2 
and  5.)  Endowing  each  ring  with  a  single  bending  stiffness  is  equivalent  to 
considering  the  ring  to  be  without  torsional  stiffness  and  attached  to  the  shell 
by  a  frictionless  bond  which  maintains  contact  but  allows  the  shell  to  slide 
freely  under  the  ring.  The  equation  of  equilibrium  may  be  solved  by  the 
Galerkln  method  in  a  manner  similar  to  that  of  references  7  and  8. 

The  simple  support  boundary  conditions  considered  herein  are  satisfied  if 
the  solution  is  taken  as  the  Infinite  series  deflection  function,  that  is, 

00 

w  =  sin  ^  ^  ^ 
m=l 


where  the  coefficients  am  are  to  be  determined.  When  equation  (A2)  is  sub¬ 
stituted  into  the  equilibrivun  equation,  the  Galerkin  method  of  solution  yields 
the  equation 


^ap 


+  M  H  y  sin  y  aj^  sin  = 

L  [_,  n  +  lZ_.^  n  +  1 


i=l 


k=l 


vhere 

p  =  1, 2, 5, 


%=r 


MxMy 


1 _ 2p^2  ^  _%P - ^ 


_ p^Z^ _ 


+  kxp‘ 


(A3) 
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d 

ring  spacing 

n 

number  of  rings 

H  = 

Dx 

P  = 

l/a 

kx  = 

.  NxL2 

n^Dx 

z  = 


R  m 


7  = 


Erir 

Dyd 


Employing  the  identity 


n 


V  • 

)  sin  — ^ 
/  .  n  + 


pjti  .  kjti  n  +  1  ^ 

-  sin  - -  =  — ~ 

1  n  +  1  2 


i=l 


vhere^  if  p  -  k  is  a  multiple  of  2(n  +  1)^ 

^npk  =  +1 

if  p  +  k  is  a  multiple  of  2(n  +  l), 

^npk  =  -1 

and  if  neither  or  both  are  true^ 

^npk  =  0 

equation  (A3)  becomes 


V 

MpUp  ^  ^  ^k^npk  “  ^ 

k=l 


where  p  =  1^  2_,  3^  ....  Equation  (Ak)  gives  the  criteria  for  buckling  if 
the  determinant  of  the  coefficients  vanishes^  that  is^ 

Aij  I  =  0  (A5) 
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The  determinant  obtained  by  expanding  equation  {Ak)  can  be  factored  into  minor 
determinants  which  correspond  to  the  buckling  modes  of  general  and  panel  insta¬ 
bility.  These  expressions  are: 


(l)  General  instability  for  a  finite  number  of  rings: 

|Aij|  = 


Mp6ij  +  (-1)  H 


=  0 


(a6) 


where ^  when  i  is  odd 


p  =  (i  -  l)(n  +  l)  +  s 


and^  when  i  is  even 


p  =  i(n  +  l)  -  s 

and  the  minor  determinant  niimber  s  equals  1,  2,  3,  .  .  .  n. 

(2)  General  instability  for  an  infinite  number  of  rings  which  corresponds 
to  orthotropic  theory  where  the  ring  stiffness  is  added  to  the  stiffness  of  the 
cylinder  wall: 

Mp  +  H  =  0  (at) 


where  p  =  1^  2,  3^  .  .  .  . 
(5)  Panel  instability: 


Mp  =  0  (a8) 

where  p  =  j(n  +  1)  and  j  =  1,  2,  3^  •  •  •  • 

The  calculations  required  to  solve  the  instability  expressions,  equa¬ 
tions  (a6),  (at),  and  (A8),  for  the  critical-stress  coefficient  kx  were  made 
for  a  large  range  of  orthotropic  cylinder  stiffnesses.  For  the  calculations 
presented  herein  a  typical  set  of  cylinder  stiffnesses  was  chosen.  These  values 
are  given  in  the  figures  and  tables  at  the  back  of  the  paper.  In  order  to  cal¬ 
culate  the  critical- stress  coefficient  associated  with  each  instability  mode, 
the  stability  expressions  are  minimized  with  respect  to  the  parameter  p,  which 
is  a  function  of  the  cylinder  length  and  the  number  of  circiomferential  waves  q 
in  the  buckling  mode.  Theoretically  q  must  be  zero  or  an  integer  number 
larger  than  one.  However,  negligible  errors  are  usually  involved  when  the  sta¬ 
bility  expressions  are  minimized  with  respect  to  p  instead  of  with  respect 
to  q,  unless  q  is  a  small  number.  The  value  of  q  was  found  to  be  suffi¬ 
ciently  large  so  that  negligible  errors  resulted  in  the  calculations  presented. 
The  required  numerical  calculations  and  minimizations  were  quite  lengthy,  and, 
therefore,  were  performed  on  a  high-speed  digital  computer.  The  exact  critical- 
stress  coefficient  was  not  found,  but  increments  of  p  were  taken  suitably 
small  so  that  only  slight  errors  resulted - 
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The  res\ilts  obtained  by  solving  equation  (A6)  are  presented  in  the  n  =  1 
through  six  blocks  and  n  =  1  through  five  blocks  of  tables  I  and  II,  respec¬ 
tively.  A  representative  exanqjle  of  the  minimization  of  equation  (A6)  for 
three  rings  is  given  in  figure  4.  Examination  of  figure  4  shows  that  for  the 
curve  for  s  =  1  there  are  two  minimums  resulting  from  a  mode  change.  The 
curve  for  b  =  2  has  a  single  minimum,  while  the  curve  for  s  =  5  produces 
the  absolute  minimum  and  is,  therefore,  the  value  recorded  in  table  I.  The 
types  of  results  shown  in  figure  4  did  not  occur  for  every  case,  but  the  range 
of  p  for  the  calculations  was  taken  to  cover  a  large  enotigh  interval  so  that 
the  possibility  of  double  minimums  was  checked. 

The  axial  buckle  shape  associated  with  the  minimums  of  figure  4  are  of 
Interest,  and  are  given  in  figure  5.  The  equation  for  the  buckle  shape  as  well 
as  the  values  of  kx  and  p  associated  with  the  buckle  shape  are  given  in 
the  figure. 

Most  of  the  solutions  for  the  buckling  coefficients  presented  herein  were 
made  by  use  of  4th-order  detenninants .  A  lOth-order  determinant  was  used  in 
selected  solutions  to  check  the  convergence  of  the  calculations.  The  lOth-order 
determinant  calculations  did  not  change  the  values  of  kx  in  tables  I  and  II 
by  more  than  0.01  from  that  obtained  by  using  the  4th-order  determinant,  so  that 
the  convergence  obtained  with  a  4th-order  determinant  was  considered  good. 

The  results  obtained  by  solving  equations  (A7)  and  (a8)  are  given  in 
tables  I  and  II  under  the  column  headings  of  n  =  00  and  panel  instability, 
respectively.  Representative  examples  of  the  minimization  of  equations  (Ay) 
and  (a8)  are  given  in  figures  6  and  7,  respectively.  The  mode  shapes  involved 
are  depicted  in  the  upper  right-hand  corner  of  the  figures.  The  results  are 
typical  of  those  found  in  maiiy  shell-buckling  problems.  Equations  (A7)  and  (a8) 
are  equivalent  to  the  equation  derived  by  Stein  and  Mayers  in  reference  5  if 
the  appropriate  wall  stiffnesses  are  employed  in  the  Stein  and  Mayers  equation 
when  general  and  panel  instability  modes  of  buckling  are  computed. 
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Figure  if.-  Minimization  of  general  instability  equation  for  three  rings.  Z  *  452.8; 


Figure  6.-  Minimization  of  general  instability  equation  and  corresponding  buckling  modes 
for  an  infinite  number  of  rings.  Z  =  k^2.Q;  y  =  62.5. 
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*'The  aeronautical  and  space  activities  of  the  United  States  shall  be 
conducted  so  as  to  contribute  ...  to  the  expansion  of  human  knowl¬ 
edge  of  phenomena  in  the  atmosphere  and  space.  The  Administration 
shall  provide  for  the  widest  practicable  and  appropriate  dissemination 
of  information  concerning  its  activities  and  the  results  thereof.” 

— National  Aeronautics  and  Space  Act  of  1958 
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